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Abstract 

We report results on the temperature dependence of the susceptibilities of 
a set of MBE-grown short-period EuTe/PbTe antiferromagnetic superlattices 
having different EuTe layer thicknesses. In-plane and orthogonal susceptibil- 
ities have been measured and display a strong anisotropy at low tempera- 
ture, confirming the occurrence of a magnetic phase transition in the thicker 
samples, as seen also in neutron diffraction studies. We suggest that dipo- 
lar interactions stabilize antiferromagnetic long-range order in an otherwise 
isotropic system and we present numerical and analytical results for the low- 
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temperature orthogonal susceptibility. 
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I. INTRODUCTION 



Magnetic films and multilayers have been a subject of intense study since they provide 
experimental realizations for various two-dimensional (2D) magnetic models.0'1 Most of the 
recent literature has focused on metallic magnetic structures. However, insulating anti- 
ferromagnetic structures provide an opportunity to study magnetic long-range order in 2D 
layered systems of localized spins. Among those, the EuTe/PbTe superlattice (SL) structures 
are of special interest since only one of the two components, EuTe, is magnetic. 

Bulk EuTe is a type-II antiferromagnet of the family of Europium chalcogenides, with 
the structure of NaCl. Its magnetic moments arise from the strongly localized 4/ elec- 
trons of the Eu 2+ atoms which are in a symmetric 8 SV/2 ground state. Thus, the europium 
chalcogenides have long been considered ideal realizations of isotropic Heisenberg models.i'S 
Antiferromagnetic resonance experiments have confirmed that, once dipolar interactions are 
taken into account, the residual anisotropy in EuTe is negligible .!'! Bulk EuTe has a Neel 
temperature (T/v) of 9.8 K. Below T/v, spins belonging to a single (111) plane are parallel 
but antiparallel to spins in adjacent (111) planes. The magnetic properties of EuTe are 
described by a Heisenberg Hamiltonian with two exchange constants: J\ (nearest neighbors, 
ferromagnetic) and J2 (next-nearest neighbors, antiferromagnetic): 

"Hex — J\ /] • Sj + J2 ^ Sj • Sj . (1) 

nn nnn 

The Sj vectors denote Eu 2+ spins, which have magnitude 7/2. J\ and J2 are not known very 
precisely; currently accepted values are: J\/k B = — 0.04±0.01fT and = O.lSiO.OLfT, 

where kg is the Boltzmann constant.!'! 

Although the Hamiltonian in eq. (p]) has full rotational symmetry, neutron diffraction 
experiments show that the spins lie in (111) planes.i This easy-plane anisotropy is properly 
accounted for by adding dipole-dipole interactions to the exchange Hamiltonian in eq. ([!]). Ull 

In this article we present experimental and theoretical studies of the susceptibilities 
of short-period EuTe/PbTe SL's. We use the notation EuTe(^)/PbTe(^) to denote a SL 
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structure whose period consists of £ monolayers of EuTe and r] monolayers of PbTe. Each 
sample used in the present work consisted of 400 such periods and was prepared so that 
rj = 3£. Since samples are grown in the (111) direction, the structure within the EuTe 
monolayers is that of a triangular lattice and the monolayers are stacked according to the 
ABC sequence. Section II describes our samples and presents results for their temperature- 
dependent zero-field susceptibilities. Section III presents a discussion of the experimental 
results, then Section IV presents a mean-field determination of the order parameter, a Monte- 
Carlo simulation of the susceptibilities of an EuTe(3)/PbTe(9) sample, and a calculation of 
the out-of-plane susceptibility at low temperature. Section V summarizes our findings. 



II. EXPERIMENTAL RESULTS 
A. Sample Characterization 

Experiments were conducted on SL samples of EuTe(£)/PbTe(3£) for 1 < £ < 7. In each 
sample the SL stack was grown on a 3000 A PbTe (111) oriented buffer layer, itself grown on 
a BaF2 (111) substrate. A 500 A PbTe cap layer was used to prevent oxidation of the highly 
reactive EuTe. Details of the MBE growth process have been published elsewhere. 0Hll The 
SL's have approximately square wave composition modulation, as reflected by the multiple 
narrow SL peaks of the high-resolution x-ray diffraction data.0 Electron spin resonance 
(ESR) experiments show very little interdiffusion at the EuTe-PbTe interface-El Furthermore, 
careful in situ scanning tunneling microscopy (STM) investigations have shown that the 
PbTe and EuTe heterointerfaces are quite smooth on a length scale of at least 200 A, with 
imperfections strictly limited to single monolayer steps.0 The ex situ cross-sectional TEM 
images also exhibit smooth PbTe and EuTe interfaces.El 

The magnetic properties of the SL's used have been studied previously by SQUID 
magnetometryQ and elastic neutron scatteringJll Magnetic hysteresis curves and neutron 
diffraction spectra taken at 1.8K and 4.2K show that for all samples with ( > 3 a transition 
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to a low-temperature ordered phase takes place at Tjy > 4.2K, the order being that of a 
type II antiferro magnet, i.e. identical to that of bulk EuTe. Since the MBE samples are 
grown along the (111) direction, this implies that the spins in each EuTe monolayer order 
ferromagnetically and are antiparallel to those in neighboring monolayers. Moreover, the 
spins lie within the EuTe monolayers. Static magnetization measurements taken parallel to 
the SL plane show no detectable in-plane anisotropy.El 

B. Experimental Conditions 

Our susceptiblility measurements used a Quantum Design (MPMS5) AC susceptometer. 
Susceptibilities were measured in two geometries : Xin was measured with the AC probing 
field h parallel to the SL plane, Xout was measured with h orthogonal to the SL plane (see 
Fig. 0). 

For each sample, Xin an d Xout have been measured as a function of temperature from 15 
K in the paramagnetic region to below the transition temperature, in a nominal zero external 
static magnetic field (i.e., less than 2 Gauss) with a 20 Hz AC probing field of 4 Gauss. Since 
the probing field is small, diamagnetic contributions from the BaF 2 substrate and the PbTe 
buffer layer can safely be ignored. A study of the dependence of the susceptibility on h and 
on the AC frequency indicates that our measurements are always in the linear regime of the 
static susceptibility. 

C. Experimental Susceptibilities 

Figure || displays the temperature dependence of the in-plane and ou-of-plane magnetic 
susceptibilities Xin an d Xout 

for samples EuTe(2)/PbTe(6) through EuTe(7)/PbTe(21). The 
experimental susceptibilities have been normalized by the number of Eu 2+ spins, determined 
using high temperature susceptibility data; no demagnetizing corrections have been applied 
to the data. In the following we will be mostly interested in Xout, whose magnitude changes 
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little below the temperature at which Xout and Xin separate. Since there is no observable dif- 
ference in Xout and x%n above this temperature we can assume that demagnetizing corrections 
are not significant. 

For all samples with £ > 1 a plot of xT 2 versus T showed that, as the temperature is 
lowered, the Curie behavior of Xin and Xout persists until they separate at a temperature 
T S .E1 Fig. |] shows such a plot for sample EuTe(3)PbTe(9). Generically, below the point of 
separation a very anisotropic behavior of the susceptibility is observed, with Xout < Xin, an 
unusual feature. Furthermore, while Xin displays a peak, Xout has a very mild temperature 
dependence below T s . In addition, the magnitude of Xout at low temperature, of the order of 
1.10 -24 emu/spin, has little dependence on the thickness of the EuTe layer. Since we do not 
have other experimental data (e.g. specific heat) that would allow us to define more precisely 
the transition temperature T/v, we assumed that it lies in between T s and the temperature 
at which Xin is maximum. T/v increases with the thickness of the EuTe layer, as expected, 
and reaches values higher than that for bulk EuTe for samples with £ > 5, an unexpected 
result. A possible explanation could be that coupling constants have values different from 
the bulk ones, a point developed further in section [V]. We comment now on samples that 
show non-generic behavior. 

1. EuTe(l)/PbTe(3) 

Xin and Xout as well as XinT 2 and XoutT 2 have been measured for T > 1.8K. They 
are plotted in Fig. |3|. One sees that in this temperature regime sample EuTe(l)/PbTe(3) 
remains in a paramagnetic phase : Xin and Xout coincide and exhibit a Curie behavior. The 
absence of a transition in this temperature range can be understood by noting that in a 
monolayer geometry the J2 exchange coupling is not present and the energy scale is set by 
the nearest-neighbor exchange coupling only, which is very small. We expect however a 
transition at a lower temperature. 
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2. EuTe(2)/PbTe(6) 



Elastic neutron scattering spectra show no long-range order at 4.2K but they do show a 
peak corresponding to type II antiferromagnetic ordering at 1.8K. Thus the broad maximum 
in Xin does not signal a transition from the paramagnetic to the antiferromagnetic phase. 

3. EuTe(6)/PbTe(18) 

The sharp drop in Xout around 5K is reproducible. However we only had one 
EuTe(6)/PbTe(18) sample. EuTe(6)/PbTe(18) is the only sample in our series that dis- 
plays this feature, as yet unexplained. 

III. DISCUSSION 

A qualitatively different behavior for Xin is expected a priori, depending on whether the 
number of EuTe monolayers per SL period is odd or even. In the first case, each period 
should behave as a ferromagnet and there should be a peak in Xin whereas in the latter case 
each period should behave as an antiferromagnet and we should see a smooth maximum. 
Although we use integers £ and rj to label our samples, the average thickness of, say, the 
EuTe layer in an actual SL period as determined by X-ray scattering is fractional because 
of the interface structure so that the odd/even effect is in fact expected to be blurred.0 

In our view, the main issues raised by our data is the existence of a phase transition in 
samples with £ > 3 at temperatures comparable to the bulk T N . As mentioned above, 
neutron diffraction spectra unambiguously demonstrate that these samples have a low- 
temperature ordered phase. Also, anisotropy in the Hamiltonian for Eu 2+ spins is negligible, 
so that our samples can be considered as representing a model 2D Heisenberg system. In such 
a system, with isotropic exchange couplings, the transition should occur at T=0K.@ The 
Hamiltonian is however incomplete and we need to supplement the exchange couplings with 
dipole-dipole interactions. Although the dipole-dipole coupling may be weak, it breaks the 
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rotational symmetry and is long-range, which prevents application of the Mermin- Wagner 



they compete with uniaxial anisotropy.EiE3 

Early workll'il has demonstrated the possibility of a phase transition driven by dipolar 
interactions in a 2D isotropic spin system. More recent work has investigated the dependence 
of Tjy on the magnitude of the dipolar coupling.i!~i3. It is our hypothesis that dipolar 
interactions are responsible for stabilizing magnetic long-range order in EuTe/PbTe SL's. 
In the next section we discuss some implications of this hypothesis. 

The observation that Xout < Xin f° r ah samples can be understood by noting that there 
are actually three pertinent susceptibilities in our system : along the direction of the 
order parameter in the SL plane, x±i i n the SL plane, but orthogonal to the direction of 
the order parameter, and Xout, orthogonal to the SL plane. In our case, x± an d Xout will be 
different because of dipolar interactions, and Xout will be smaller than x±- Usually one has 
X\\ — X-L) an d in our case we expect x\\ < Xout < X±- The assumption that in each atomic 
layer the spins belong to domains with random orientations yields : Xin — iOftl + X±)- It is 
thus possible for Xout to be smaller than Xin- 

IV. THEORY 

In this section, we use the following Hamiltonian to describe the interaction of Eu 2+ 
spins: 



where the third sum runs over all sites i and j in the SL and on spin components a, (3. Here 
Qap(?j — is the dipolar tensor which reads: 



In this expression, /i& is the Bohr magneton and g is the Lande g factor which we will take 
equal to 2. Throughout this section we will approximate the spins Sj by classical vectors. 




TC — J\ Sj • Sj + J 2 Sj • Sj + ^2 Qa/3{ r j ~ r i)S"Sj , 



(2) 





(3) 
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This is justified by the large magnitude S = | of the Eu 2+ spins and the fact that we don't 
expect quantum effects in the temperature range we will be considering. The structure 
within the EuTe layers is that of a triangular lattice and the layers are stacked according to 
the ABC sequence. We will take the x and y axes in the plane of the layer and the z axis 
orthogonal to the layer plane. 

A. Orders of magnitude 

Let us denote by E\\ and E± the dipolar energies per spin of an EuTe monolayer, assuming 
the spins are ferromagnetically aligned either in the layer plane or orthogonal to the layer 
plane. E± is given by : 

*i =£««('<) = ^Ei, (4) 

i A i 'i 

where the sums run on the sites of a single layer. Eh is given by : 

i ^ i ' i 

where we have used the fact that Q is traceless and rotationally invariant in the xy plane. 
Notice that En < 0, which favors in-plane alignment, as observed in neutron diffraction 
experiments. Using the value of J^i ~ir — 11.035/a 3 for a triangular lattice and the value of 

z 

the in-plane lattice parameter a = 4.6A taken from X-ray data, we get En = — 0.86K and 
E± = 1.72K, to be compared with the exchange energies J\S 2 = — 0.5K, J 2 S 2 = 1.8K and 
the average exchange energy per spin E ex = —6J 2 S 2 = 11K, which we have estimated using 
bulk EuTe values for J%,J 2 . 

Next we estimate the size of the interlayer dipolar energies. To this end we have used 
Ewald summation techniques that allow one to rewrite \ sums as fast-converging series.ilil 
Let us consider two neighboring EuTe monolayers (1) and (2) a distance h apart. Let r 
be the vector joining a lattice site in layer (1) to a lattice site in layer (2). Assuming that 
the spins in layer (1) and (2) are all ferromagnetically aligned but with opposite directions 
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depending on which layer they belong to, the interaction energy of a spin in layer (1) with 
all spins in layer (2) is: 

S=-^^E§^«(G') (6) 

where the sum runs over all reciprocal lattice vectors G, and G and G x are, respectively, 
the modulus and the x-component of G and A is the area of the triangular lattice unit cell 
in the layer. We thus get : 

£ = 0.1014 i^bSf = Q317K 

This energy is much smaller than the intralayer energy E\\. If the spin density in the layer 
were uniform, no field would be created outside the layer and this energy would be 0. The 
very existence of a lattice structure within the layer makes it finite. Furthermore, the 
interaction energy is expected to decay fast as the distance from the layer becomes larger 
than the in-plane lattice constant. For instance, at a distance two layers away, this energy 
is —8.8 10~ 5 K, and three layers away it is 4.7 10~ 7 K. 

As a result, we can safely discard all interlayer couplings as well as interperiod couplings, 
which couple spins belonging to different SL periods. Although small, the intralayer coupling 



has to be retained for the reasons mentioned in section III 



B. Mean Field analysis 

We now turn to a Mean Field treatment of our problem. Our aim here is to identify 
the order parameter for the phase transition rather than find the expression for T/y.. We 
consider a single period of a SL which consists of N EuTe monolayers. The spins are labelled 
with two indices: i denotes their position in a layer, n the layer to which they belong. We 
rewrite the Hamiltonian as: 

7~C = J]] Hi na j m [3 S in Sj m (7) 

ijnm,cef3 

where the matrix H is defined by: 
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H ina>jm p = 8 a ? Jfo -r h m-n) + 5 nm Q a(3 (r 3 - r<) (8) 

in which J denotes the matrix of the exchange couplings and its elements are equal to J\ if 
(in) and (jm) are nearest neighbors, to J2 if they are next-nearest neighbors, and to zero 
otherwise. The S nm factor expresses the fact that interlayer dipolar couplings are neglected. 

A mean-field calculation in our context amounts to diagonalizing H and finding its lowest 
eigenvalue. The magnitude of the latter determines T/v while the associated eigenvector 
defines the order parameter for the transition.il 

Since Q is diagonal in the layer indices we concentrate first on J. We define in-plane 
Fourier transforms S qn for the spins through: 

Sqn = — jjj E Si„e jq ' r * (9) 



where M is the number of spins per layer. Likewise we define J(q, m — n) as: 

J(q,m-n) = J2 </(d, m - n)e iq d . (10) 

d 

Since we expect the ordered phase to be homogeneous in the plane of the layers, we now 
restrict ourselves to the q = sector of the Hamiltonian. The neighbors and nearest- 
neighbors of a spin in layer n all belong to layers n — 1, n, n + 1. The only non-zero matrix 
elements of J(0, m — n) are thus: 

J(0,0) = 6Ji and J(0, 1) = J(0, -1) = 3(Ji + J 2 ). 

The q = part of the Hamiltonian now reads: 

7"^q=0 = E A nm SonSom (H) 

nm 

where A nm = J(0,m — n). Now A can be diagonalized in the basis of N- dimensional 
orthonormal vectors Tk, whose components are: 



where k is an integer ranging from 1 to N. The corresponding eigenvalues are: 
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a k = 6Ji + 6(Ji + J 2 )cos [jj—j J . (13) 

In the same way as for J in eq. ([0]) one can define a Fourier transform Q a p{<\). It is 
diagonal for q = with: 

Q XX (0) = Qyy(0) = -C (14) 

Qzz(0) = 2C (15) 

where 

c = ^E^ (is) 

i ' i 

The lowest eigenvalue of matrix H is thus a^r — 2C and the mean field T N is given by: 

2S" 2 

T JV = --— (a^-2C) 
3 k B 

= S(- 6Jl - 6(Jl + J2)cos (/^Tl) + 2C ) ' (17) 

Using bulk values for Ji, J 2 we find T/v = 6.93K. The associated order parameter is a linear 



combination of the in-plane projections of the spins with weights {Tk)n defined in (0): 

m,^x>G^) 5 „ m 

where S n denotes the in-plane projection of the sum of all spins belonging to layer n. Note 
that, as N — > oo the usual antiferromagnetic staggered magnetization is recovered. We will 
in the next section check that is indeed the correct order parameter. 



C. Monte Carlo simulation 

We have performed a Monte Carlo simulation for our system in order to check that dipolar 
interactions can drive a transition at a temperature T/v > and can generate an anisotropy 
in the susceptibility similar to that observed in experiments. We have also checked the 
relevance of the mean field order parameter. 
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We present here results of a Monte Carlo simulation carried out on a system of 3 layers, 
each consisting of 23x23 spins, with periodic boundary conditions in the plane of the layers, 
in order not to introduce in-plane anisotropy. Each layer is thus mapped to a torus. The 
full Hamiltonian (|7]) has been used, where the distance between pairs of sites has been 
taken to be the smallest distance on the torus between sites i and j. All couplings in the 
Hamiltonian have been expressed in units of J\ and we have used EuTe bulk values for 
the ratios Jij J\ = —3.75 and (g/if,) 2 /( Jia 3 ) = 0.64, which are the only parameters of our 
model. The heat-bath algorithm with sequential updating of the spins has been used, with 
400 equilibration sweeps and 2000 sweeps with a measurement after each sweep. Error bars 
have been carefully computed as standard deviations of estimators of the observables. Three 
susceptibilities have been evaluated: Xz, along the normal to the layers, corresponding to 
Xout, an d Xx an d Xy m the plane of the layers. The specific heat has also been measured 
to check that it has a limit of ks per spin as T — > 0, a general property of classical spin 
systems.0 It is shown in Fig. ^ for the EuTe(3)/PbTe(9) system. Susceptibilities have been 
computed as: 



where S is the total spin and J\f is the number of spins per layer. 

Because the simulation has been done on a finite system, the in-plane rotational symme- 
try is not broken in the ordered phase with the result that the above averages are not well 
defined. We have thus evaluated Xx and Xy only at temperatures higher than the transition 
temperature. 

Experimental and simulated susceptibilities have been plotted as a function of temper- 
ature in Fig. p. One can see that the simulation qualitatively reproduces the anisotropy 
observed in the experiments, with a remarkable flatness of Xz (Xout) at low temperature. We 
have also run simulations on 3-layer systems of sizes 29x29 and 13 x 13 and these two systems 
didn't show any significant difference in the transition temperature or the low-temperature 
magnitude of Xz- 
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We then checked the relevance of the order parameter M^ found in the mean field 
approach of section B. In the same way as Mjy in eq. flIB|) has been defined using vector 
Tk=N each vector Tk can be used to build an order parameter M&, defined by: 



M* 



^vi^hh (19) 

using the same notations as before. The M^'s are linear combinations of the in-plane 
projection of the spins. One could similarly define linear combinations of their z components, 
although we know that spins order in-plane. For a 3-layer system, three order parameters 
can be defined: 



- - 



Mi = — — Si + S 2 + — <S 3 



M 2 = i= (5! - S 3 ) 

One can define susceptibilities Xk for the moduli of these three order parameters k = 1,2,3 



as: 



V/- T77^ (<M£> - (|M fc |) ; 



NMk B T 

Because vectors Tk are normalized, the Xk ai l have the same leading behavior at high tem- 
perature. These susceptibilities are plotted in FigM X3 is the susceptibility that displays a 
sharp peak, thus suggesting that M3 is the appropriate order parameter for describing the 
phase transition. 



D. Low-temperature orthogonal susceptibility 

Since Xout is flat at low temperature, it is desirable to have an estimate of its magnitude. 
We present here a calculation of Xout using the spin Hamiltonian given in eq. (0). 

T~t = Hex + 7~Ld 

where the exchange term is: 
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Hex = J ( V 3 ~ r *> m ~ n ) ' S i™ ( 20 ) 

ij,nm 

and the dipolar term is: 

H^^Q^vj-v^S^. (21) 



ij,n 



In each layer we define a frame of reference as shown in Fig. |8|, such that the y axis lies 
along the direction of the layer magnetization, while the direction of the z axis, orthogonal 
to the layers, is the same for all layers. The spin components af n in the layer- dependent 
frames are related to the Sf n through: 



qx 
^in 


= (-1)X„ 




qy 

'-'in 


= (-l)Vn 


(22) 


qz 

^in 


= °in ■ 





At low temperature a\ n and af n will be small, while <rf n will be finite, with a fixed sign. Xout 
is defined as: 

Xout ~ NAfk B T {[ ' ' 
where a z = 2i n of n . F° r classical 3D spins ((cr z ) 2 ) is given by: 



<(0 2 > = 7/n^f ipfe-P" (23) 

Z J in S \ a in\ 

where (3 = The partition function Z reads: 

j y s Ian 



At low temperature (cr.f n ) 2 and (crf n ) 2 will be of order k B T if \out is to be finite. We can thus 
let the integrals run from — oo to +oo, instead of — S to +S. We note that a y in appears both 
in the integration measure and Ti and is given by: 



< = ^S 2 -(af n 2 + af n 2 ). 
We now expand a\ n in powers of [p\ 2 + crf n 2 ) keeping only the lowest order terms. This 



amounts to an expansion in powers of the temperature. The measure becomes da z n daf n /S 
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In the expression of /37i only the first-order term need be kept. We are then left with a 
quadratic form for af n and af n which we need to diagonalize in order to calculate ((a z ) 2 ). 
Let us first consider TC ex - In our approximation Sj n • Sj m reads: 

/ z2 + X 2 z 2 + x 2 X 

Sin ■ S jm = a\y jm + (-ir~ n [oZofn + (S- a ™ ){S - Jm 2S 3m ) j • (24) 

After discarding constant terms, we get the following expression for 7i ex 

K ex = E ( J (*3 -n,m-n) + aJijSnJ (a z m a z jm + (-l) m - n af n ^ m ) (25) 

ij,nm 

where a n = 3(J2 — Ji) if n = 1 or n = N and a n = 6J2 otherwise. Let us define Fourier 
transforms for of n , <rf n and J(r^ — r^m — n) in the same way as in @ and (|T0|). We now 
have: 

W« = E (^(q,m - n) + oO™) (< n *< m + (-1) B <C(-1) W <C) ■ (26) 

mn,q 

Let us now turn to 7^^ and expand of n . Because of the layer geometry Q xz = Q yz = 0. 
Furthermore, since Q(q = 0) is diagonal, non-diagonal terms Qy X (*j ~ r i) <T in (T jn d-O not 
contribute at the quadratic order, but rather yield a term linear in af n which vanishes when 
summed on j. We are thus left with diagonal terms only. The contributions of Q xx and Q zz 
are then: 



E^(q)|<J 2 + Q ra (q)l<J 2 (27) 

qn 

and that of Q yy , after expanding af n to first order, is: 

-Q ra (q = o)E« 2 + < 2 ) = cE(l<nl 2 + l<J 2 ) 



qn 



Since only the q = mode contributes to a z , we now restrict Tid to its q = part. After 
using ( |L4D and (|1|) we get 



^ = 3CE^on 2 (28) 

n 

and 
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H ex = £ (J(0,m - n) + a n 5 nm ) (a z 0n a z 0m + (-l) B ^(-l) m a£J . (29) 

ran 

We only need diagonalize the z part of Ti. which reads: 

Hz = B nm ffon ff Om ( 30 ) 

nm 

with 

B nm = (3C + a n )8 nm + J(0,m-n) (31) 



Matrix B has the same form as matrix A used in section [IV B| with the difference that the 

matrix elements at both ends of the principal diagonal are different. B can be diagonalized 

in the basis of orthonormal vectors U k defined by: 

1 / 1 kir\ 
(t4)n = ^sin(^--)-j (32) 

where the integer k ranges from 1 to N. The normalization factor is equal to y/~N if 



k = N and \ N/2 otherwise. The corresponding eigenvalues are: 



h = 3C + 6{Ji + J 2 ) f 1 + cos(— ) J . (33) 



Expanding a z on the basis of the a^ k defined by: 

a ok = ( U k)n cr z 0n 

n 

we get: 



where 



if A; even 
Afc = { _ (34) 



^ if k odd. 



The integration over the of n in eq. (^3|) is straightforward since the a^ k have Gaussian 
weights and we get the following result for Xout at low temperature: 

X ° Ut 2N* t 3C + 6( J x + J 2 ) (l + cos(^ m 
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where the weights fk are given by: 

if k even 

fk - 



2 

sin^ff) 



if k odd and k ^ N (36) 



. if k odd and k = N. 



The energy C is defined in ([16]) and equals 2.759 . The /& satisfy the following sum 
rule: E k fk = N 2 . 

The Xout we have found is temperature- independent; it is in fact the first term in an 
expansion of Xout in powers of T. Using bulk values for the couplings, we find that for a 
3-layer system \ out = 0.025 in excellent agreement with the simulation result in Fig. 
|6]: (0.026 ± 8.10~ 4 )^4jt-. As the number of layers is increased, Xout slowly decreases. For a 



7-layer system, expression fl35|) yields: Xout = 0.017 



IJil 



V. DISCUSSION AND CONCLUSION 



The simulation results reported in section |1V C| for a 3-layer system qualitatively re- 



produce the anisotropy in susceptibilities observed in the experiments. However they do 
not agree quantitatively with the experimental results. We list below the experimental and 
theoretical values of T N and Xout at low temperature, in dimensionless units. 

T^ xp = 16.1 T™ c = 9.5 ± 0.5 in units of 

X e Z = 0.013 xZt = 0-026 ± 8.10- 4 in units of « 
Our Monte Carlo simulations have been run on finite systems: although Tjv and Xout didn't 
change appreciably when we increased the size of the system to 29x29 or decreased it to 
13x13, we cannot rule out finite-size corrections. For Xout, however, analytic and simulation 
results are in excellent agreement which indicates that finite-size corrections are not signif- 
icant in the low-temperature region. The numerical values listed above suggest that the 
effective value of J\ in the samples is approximately twice as big as in the bulk. However, 
since the simulation uses bulk values for the ratios Jij J\ and {g^) 2 / [a? J{) this would imply 



that J2 and the dipolar coupling are rescaled by the same factor. This cannot be the case, 
as the value of the dipolar coupling only depends on the in-plane lattice parameter, known 
from X-ray spectra. One possible explanation is that the exchange constants are different 
in the SL's relative to bulk values. Because of a 2.1% lattice mismatch between EuTe and 
PbTe, the SL will be strained and the in-plane and out-of-plane lattice constants will be 
different from one another, and different from their bulk values^. As a result, the exchange 
couplings will also be somewhat different. Within the family of Eu chalcogenides the lattice 
constant increases as the size of the anion increases from O to Te which makes it possible 
to study the dependence of the exchange couplings on the lattice parameter.0 In our case 
the in-plane lattice constant is reduced with respect to the bulk value while the out-of-plane 
constant is increased. Thus, in the SL's, the in-plane J\ is likely to increase, while the 
out-of-plane J\ will decrease. 

Although dipolar interactions account for the flatness of Xout at low temperature, we can 
expect a single-ion anisotropy term of the form K^ in Sl^ to have the same effect. Let us 
replace Tid with such a term. The calculation of \out is similar, except that in equation 
3C is replaced by k. Thus for single-ion anisotropy to have the same effect as dipole 
interactions, k would have to be of the order of 0.2K, which is larger than the exchange 
couplings, an unlikely situation. 

In conclusion, we have presented experimental susceptibilities of EuTe/PbTe short-period 
antiferromagnetic superlattices. We suggest that dipolar interactions may stabilize long- 
range order in these 2D structures. Additional theoretical work along with more precise 
susceptibility data and specific heat measurements are needed to confirm this hypothesis 
and to study the critical behavior, an aspect not touched upon in the present work. 
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FIGURES 
(a) (b) 




Xin Xout 

FIG. 1. The two principal orientations of the sample with respect to the AC probing field for 

the Xin & n d Xout (b) susceptibility measurements. 
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FIG. 2. Xm (in-plane, •) and \out (orthogonal, o) susceptibilities in emu normalized per Eu 
atom, for samples EuTe(2)/PbTe(6) through EuTe(7)PbTe(21). 



24 



5e-24 




1e-22 



10 15 20 25 
T(K) 



30 



FIG. 3. (a) Xin (in-plane, •) and \out (orthogonal, o) susceptibilities in emu normalized per Eu 
atom, for sample EuTe(l)PbTe(3). (b) Same susceptibilities, multiplied by T 2 . 
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FIG. 4. A plot of X T 2 for sample EuTe(3)/PbTe(9) in the paramagnetic phase (squares) and 
in the ordered phase, in the in-plane (•) and out-of-plane (o) directions, in emu normalized per Eu 
atom. 
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FIG. 5. Simulated specific heat for a EuTe(3)/PbTe(9) system in units of ks- 
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FIG. 6. (a) Measured Xin an d Xout for the EuTe(3)/PbTe(9) sample, (b) Simulated x,y,z 
susceptibilities for a 3-layer system of size 23x23. 
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FIG. 7. Simulated susceptibilities corresponding to order parameters for a 3-layer system 
of size 23x23. 
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